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Математическая модель динамики вращения ветротурбины и поворота малой 
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автоматической системы ориентации 
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Представлено математическое моделирование динамики вращения ветротурбины и поворота ветроэнергети- 
ческой установки. Показано распределение углов атаки лопасти по азимуту вращения ветротурбины при 
угле косой обдувки 30°. Видно, что обтекание при наличии ветра несоосного с осью вращения приводит к 
существенной неравномерности углов атаки и, соответственно, к изменению сил на лопастях в процессе 
вращения. 
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Введение. Эффективность горизонтально-осевой ВЭУ зависит от многих факторов. Одним из 
важных показателей является способность установки быстро перестраиваться навстречу ветру с 
наименьшим углом к нему. В то же время, вследствие инерционных сил, которые затрудняют пе- 
рестройку ветроустановки вслед за изменениями метеорологической обстановки, ВЭУ постоянно 
находится в состоянии присутствия косой обдувки, т. е. наличия угла между осью вращения и на- 
правлением ветра. Свою лепту в наращивание величины вносит и трение в опоре горизонтально- 
осевой ВЭУ. Для того чтобы оценить влияние различных факторов и определить эффективность 
ВЭУ уже на стадии разработки, целесообразно иметь математическую модель ветроустановки. 
Для её создания необходимо рассмотреть основные силы, действующие на ветротурбину. 
Моделирование динамики поворота ВЭУ. Как известно, вращение ветротурбины обусловле- 
но наличием аэродинамических сил, создающих крутящий момент на валу генератора. Однако 
этому вращению препятствует момент генератора, появляющийся вследствие наличия сил трения 
в мультипликаторе, опорах вала и сил сопротивления вследствие электромагнитных сил. 

Суммарное уравнение моментов на валу ветрогенератора можно представить следующим 
образом: 


04 
а = (Ты +), (1) 


где М, — крутящий момент, обусловленный наличием аэродинамических сил; М, — момент со- 
противления на генераторе, обусловленный наличием электрической нагрузки; М, “М0: 
мент, обусловленный силами трения в мультипликаторе и опорах вала; /„„ — момент инер- 
ции ветротурбины; л — передаточное отношение мультипликатора генератора; /„ — момент 
инерции ротора генератора; и, — угловая скорость вращения ветротурбины. 


При этом разворачивающий момент ВЭУ равен: 


О 

= 2 
5 ы а: -М,, — 1 СЕ / (2) 
где М, — момент от аэродинамических сил на ветротурбине, обусловленный КОСОЙ обдувкой; 
М — момент от сил, действующих на оперение, М ВЕ гироскопический момент ветротур- 
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бины; М, — момент от сил трения в опоре ветротурбины; Г, ° — момент инерции ветротур- 


2 вэу 


бины относительно оси поворота; «, — угловая скорость поворота ВЭУ относительно верти- 


кальной оси. 
Момент генератора является функцией электрической нагрузки М; =Ё (^) ‚ трение в опо- 


ре величина постоянная 


2 п 


а 
М. =К > Тв» 9 , (3) 


где К, — коэффициент трения в подшипнике; /7,., — масса поворачивающейся части ВЭУ; @, — 


ВЭУ 
диаметр сепаратора опорного подшипника. 
Моделирование динамики вращения горизонтально-осевой ветротурбины. Крутящий 
момент, обусловленный наличием аэродинамических сил, зависит от характера обтекания ветро- 
турбины и, в первую очередь, наличием косой обдувки. Он может быть рассчитан несколькими 
методами. Наиболее оптимальным, с точки зрения сложности и точности, является метод расчёта 
единичного элемента лопасти [1, 2, 3]. Этот метод расчёта основывается на допущении, что каж- 
дый элемент лопасти ветротурбины можно рассматривать как отрезок профиля, двигающийся по 
винтовой линии. Аэродинамические силы рассчитываются, исходя из результирующей скорости 
потока, обтекающего этот профиль. Причём считается, что прилегающие части лопасти не влия- 
ют на характеристики рассматриваемого отрезка профиля. Силы и крутящий момент винта полу- 
чаются интегрированием элементарных сил и моментов отдельных элементов лопасти. Соответст- 
венно задача определения аэродинамических характеристик ветротурбины решается при сле- 

дующих допущениях: 

— лопасть является абсолютно жёсткой на изгиб и кручение; 

— аэродинамические характеристики профиля сечения лопасти не зависят от скольжения 
потока вдоль лопасти и влияния центробежных сил на пограничный слой. 

На элементе лопасти при вращении возникают аэродинамические силы, обусловленные 
наличием как окружной, так и набегающего потока через ветротурбину. Представленное на рис. 1 
соотношение скоростей сечения лопасти соответствует лопасти, находящейся в верхнем положе- 
нии с обдувкой набегающего потока под углом а, к оси вращения. В этом случае вертикальная 
составляющая скорости потока И, отрицательна и угол атаки а, меньше угла установки. 

Скорость продольного движения представляет собой сумму скорости набегающего потока 
и окружной скорости. 

И’ = ог +И/ Уп, (4) 
где г — текущее значение радиуса; « — угловая скорость несущего винта; И/ — скорость ветра; 
Ц — азимут лопасти. 


Вертикальная составляющая в плоскости сечения профиля скорости элемента лопасти бу- 
дет равна: 


И, =И/ с05а, , (5) 


где а, — угол между осью вращения и направлением ветра. 
Угол притекания имеет вид: 


В=фча, (6) 
г { И, 

р (7) 
а, =В-ф (8) 
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Рис. 1. Распределение аэродинамических скоростей и углов на сечении лопасти 
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Рис. 2. Распределение аэродинамических сил на сечении лопасти 
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Значение осевой силы элемента лопасти равно: 


ат =аХ „тр + аУ„, со5В. (9) 
Значение элемента силы, направленной на вращение ветротурбины: 

АО = ау „ эптВ - аХ „„ со5В (10) 
АМ, = О (11) 

2 (/? 2 
аХ = вс Бар = ыы И (12) 

Е 2 

2 2 2 

м РЕ с Ба = НЫ Бар , (13) 
р 2 ) 


где с, — коэффициент подъёмной силы профиля на участке элемента лопасти; с, — коэффици- 


ент силы сопротивления профиля на участке элемента лопасти; И. — суммарная скорость об- 
текания элемента лопасти; И„, И, — проекции скорости движения элемента лопасти; В — 


хорда лопасти; р — плотность воздуха; (г — длина элемента лопасти вдоль радиуса. 
Величина коэффициента подъёмной силы в линейной зоне может быть определена по 
формуле: 
О, 
с, =а а, =а, [+ (14) 
Хх 


где а, — производная коэффициента подъёмной силы по углу атаки; а, — текущий угол атаки 


сечения относительно точки нулевой подъёмной силы профиля; ф — угол установки профиля. 
Текущий угол установки профиля представляет собой: 
(-Ю) 
К-К 
Принимая во внимание тот факт, что помимо вращения ветротурбины, присутствует и 
второе вращательное движение, а именно вращение всей установки ветрогенератора, тогда ско- 
рость продольного движения, а также вертикальная составляющая скорости элемента лопасти 
принимают несколько иной вид. 


Плечо скорости вращения И. , действующей на произвольную точку лопасти, равно: 


го = Г? +? , (16) 


где г, — плечо вращающей силы, Г — расстояние от центра лопасти до точки приложения сил, 


ф; =Ф, +ДФ, =Ф, +Дф (15) 


т — расстояние от оси вращения ветрогенератора до его ветротурбины. 


у = агс9- (17) 

Ир =©,Г„ С05У (18) 

И’ = г -И/ чпа, со$ф (19) 
И, =И/ со5а, +И.„, эф (20) 


Для проведения интегрирования величины аэродинамических сил ветротурбины по радиу- 
су необходимо учесть потерю эффективности ветротурбины на краю диска вследствие перетека- 
ния воздуха на концах лопастей. Вследствие разницы давления воздух перетекает с нижней по- 
верхности на верхнюю, уменьшается разница в давлении и уменьшаются и силы в области концов 
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лопастей. Зона распространения концевых потерь зависит от удельной нагрузки на лопасть и на 
всю ветротурбину. 

Однако мало нагруженные ветротурбины можно рассчитывать при принятии некоторых 
простых допущений, касающихся концевых потерь. Одним из них является допущение о том, что 
часть лопасти длиной, равной половине концевой хорды, не создаёт аэродинамических сил. В 
этом случае относительный радиус рабочей части лопасти равен: 


Ь 
В =1-—. 21 
т. (21) 
Осевое усилие на ветротурбине равно: 
360 ВК 
Т=2 | | ат. (22) 
0 Г 
Интегрированное значение момента лопасти определяется так: 
360 ВК 
М, =2 | | г,аТ, (23) 
0 [ 


где 7 — количество лопастей. 
Таким образом, момент ветротурбины представляет собои: 


360 — ВК р(( +И/ чпа, с0$ 4) + (И/ соза, + ог, со$ у 5т +] | 


М = с Баг |$тВ + 
оу 5 Е В 
| | (24) 
р((г +И/ па, со$ф)` + (И/ соза, + ®г,, со5у пу) 
С, [6088 1,67 
Крутящий момент определяется следующим образом: 
360 ВК 
М, =2 | у | Гао. (25) 
0 [ 
При численном интегрировании лопасть разбивается на элементы длиной: 
БК —Ю 
аг„ = рывкь | (26) 
п, 
где п» — количество отрезков, на которые разбивается лопасть. 
Радиус #го элемента лопасти равен: 
= +5, (27) 
ие = Гир) + [ (28) 


Выводы. На основании предложенной выше методики было выполнено моделирование экспери- 
ментальной ветротурбины диаметром 2 м на разных углах косой обдувки. На рис. 3 показано рас- 
пределение углов атаки лопасти по азимуту вращения при угле косой обдувки 30°. Видно, что 
обтекание при наличии ветра, несоосного с осью вращения, приводит к существенной неравно- 
мерности углов атаки и, соответственно, распределения сил на лопастях в процессе вращения. 
Причём если угол между направлением вращения и направлением ветра присутствует в горизон- 
тальной плоскости, то основной перекос аэродинамических усилий присутствует в вертикальной 
плоскости. 
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Рис. 3. Распределение аэродинамических углов атаки профиля ветротурбины диаметром 2 м (ветер 6 м/с, угол косой об- 
дувки 30°) 





-20,00] 
Рис. 4. Мощность на валу ВЭУ в зависимости от угла косой обдувки 
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Рис. 5. Значения моментов на ВЭУ от аэродинамических сил и сил трения 


Расчёт крутящего момента при этом показал, что косая обдувка приводит к существенно- 
му снижению эффективности ВЭУ. Так, при угле 15° снижение вырабатываемой мощности состав- 
ляет 11 %, а при 30° эта величина вырастает до 35 % (рис. 4). 

Выполненные далее расчёты поворота всего ветрогенератора при постоянной частоте 
вращения позволили определить моменты аэродинамических сил на ветроколесе и стабилизаторе 
(рис. 5). Значение момента трения выглядит как две горизонтальные линии. Это фактически зона 
нечувствительности ВЭУ к повороту. Видно, что до значения угла косой обдувки 14° ВЭУ не будет 
совершать поворот относительно своей вертикальной оси. 

Таким образом, можно констатировать, что ВЭУ всегда находится в зоне косой обдувки не 
только вследствие инерционных сил, но сил трения. 
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5УЗТЕМ 


А. 5. Синпом, У. 1. Саропоу 
(Роп Чае Тесптса! Упмег®Ку) 


Тре та тетайса!/ зити!айоп оЁ иипа игбте гоаНоп апа допао/а &игп/пд 15 ргезетеа. Тре Басе апд/е-оРаваск 
Острийоп т айтиё о! те илпа фигбте гоаНоп ии 30° апое о! [ейега! илпа 15 5роип. Тре иду $Вои/5 па 
те попаха! гоавоп Пои/ саибез сИйса! птедиагйу оЁ {Ве апд/е$ о! авасК апа арргорпаёе/ ГЮюгсе иапайоп оп те 
Ы/аде$ аитпа {Ве гоаноп. 

Кеуигога$: аегодепегаог {иг адупат/с ути!аНоп, илпа Фиготе гоайоп аупат/с эпти!айоп, Басе. 
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